NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



s£S Axiomatix 


(NAS A-CH- 160557) Ku-BAUD RADAR THRESHOLD N80-19322 \ 
ANALYSIS (Axiomatic, Los Angeles, Calif.) I 
37 p HC A03/ME A01 CSC! 171 J 

Unclas I 
G3/32 47544 





KU-BAND RADAR THRESHOLD ANALYSIS 


Contract NAS 9-15795 
(Technical Monitor: E. B. Walters) 


Prepared for 

NASA Lyndon B. Johnson Space Center 
Houston, Texas 77058 


Prepared by 


C. L. Weber 
A. Polydoros- 
Axiomatix. 

9841 Airport Blvd., Suite 912 
Los Angeles, California 90045 


Axiomatix Report No. R8001-2 
January 11, 1979 


GLOSSARY OF SYMBOLS 


1 = RF frequency index; let! • I max 3 
L = Range gate Index i {;{! 

k * Time index for a specified pair (I,L); ke[0,N-l] 

m = Doppler filter index; me[0,N-l], N = 16 

Sj(L,k), S T (L,k) = Real antf imaginary parts, respectively, of the signal 

component at the output of the Lth range gate, at 
time k. 

N R (L,k), N,(L,k) = Real and imaginary parts, respectively, of the noise 
K 1 component at the output of the Uh gate, at time k. 

PSI(L,k) = _ Real part of the total output of the Lth gate at 
S R (L,k) + N R (L,k) " time k. 

PSQ(L,k) = _ Imaginary part of the total output of the Lth gate 
Sj(L,k) + Nj(L,k) " at time k. 

F(L,m) = Output of the mth doppler filter following the Lth 
gate. 

F R (L,m), Fj(L,m) = Real and imaginary parts of F(L,m), respectively. 

2 

cr = Variance of the in-phase and quadrature phase 
s Gaussian signal components. 

2 

a = Variance of the in-phase and quadrature phase 
n Gaussian noise components. 

N = Order of the Discrete Fourier Transform (DFT) filter. 

| • | = Norm of ( 0 

p = Noise correlation coefficient. 


SUMMARY 


In this report, the statistics of the CFAR threshold for the 
Ku-band radar are determined. Exact analytical rc^ts are developed 
for both the mean and standard deviations in the designated search mode. 
The mean value is compared to the results of a simulation reported in 
[1], The analytical results are more optimistic than the simulation 
results, for which no explanation is offered. 

The normalized standard deviation is shown to be very sensi- 
tive to signal -to-noise ratio and very insensitive to the noise corre- 
lation present in the range gates of the designated search mode, The 
substantial variation in the CFAR threshold is dominant at large values 
of SNR where the normalized standard deviation is greater than 0.3. 
Whether or not this significantly affects the resulting probability of 
detection is a matter which deserves additional attention. f 

On the optimistic side, the threshold setting and target return 
are correlated; this leads us to conjecture that this variation may not 
appreciably affect the probability of detection, On the pessimistic 
side, there is a substantial variation of the CFAR threshold setting 
away from that developed from the noise-only condition. 


* f * * 


1,0 


INTRODUCTION 


The constant false alarm rate (CFAR) thresholding scheme in the 
Shuttle Ku-band radar Is analyzed for the "designated mode" of operation. 
In particular, both the mean and standard deviations are determined for 
the threshold ievel. 

In search, there are two basic modes of operation: designated 
and undesignated. In the designated mode, range being designated, there 
are two overlapping range gates of width 3t/2, where t is the transmitted 
‘pulse width. Four nonoverlapping range gates of width t are used in the 
undesignated mode. 

Sixteen pulses are transmitted at each of the five RF frequen- 
cies. When range designation is available, the pulse width and pulse 
repetition frequency are functions of the designated range. 


v 
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2.0 DESIGNATED MODE THRESHOLD 

The basic signal processing of the Shuttle Ku-band radar for 
the designated mode in search is shown in Figure 1. Only those signal 
processings pertinent to the CFAR threshold formulation are shown. For 
a more detailed description of the signal processing for the Ku-band 
radar, see [2]. 

The output of the IF filter is downconverted to a complex base- 
band waveform 


I + jQ = SI(I,J,k) + jSQ(I,J,k) 


( 1 ) 


Before A/D, the I and Q waveforms for the kth pulse are given by 


I + jQ = Aj P(t-kTp) exp|d (w d t + 0j) 


+ N c (t) + jN s (t) 


( 2 ) 


where 

A t = the random amplitude of the target return which has the 

^ 2 
Rayleigh probability density function with parameter a s , 

which represents signal power 
0j = a random phase uniformly distributed over (0,2ir) 

toj = the doppler frequency, which is neglected in this analy- 
sis. The effect of doppler on the final results is not 
expected to be appreciable. 

P(t) = pulse shape of width t seconds 
T p " (PRF)- 1 

N c (t), N g (t) = independent zero-mean Gaussian processes with one-sided 
power spectral density Ng W/Hz and one-sided noise band- 
width f c (f c = 237 kHz), which is the 3 dB bandwidth of 
Hj(f) (see Figure 1). 


The noise power in N c (t) and N s (t) is therefore given by 


N 0 f c 


C3) . 


for each process. 


Output of 



Figure 1. Signal Processing for Designated Mode in Search 
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The integration process of the presum is also shown in Figure 1 
for the designated mode in search. We assume that the received pulse is 
ideally designated so that it appears exactly between the two range gates 
of width 3t/2, as shown in Figure 1. This, coupled with neglecting the 
doppler effect, maximizes the effect of the signal received from the tar- 
get on the CFAR threshold. 

With these assumptions, the presum output for the jkth pulse 
for the Lth range gate is 


PSI(L,k) + 3 PSQ(L,k) « 


t+3r/2 




(t) + j SQ(t) 


dt 



exp 


t s +3-r/2 

We, )+/ 


N c (t) + J N s (t) dt 


( 4 ) 


where the sum of the samples at the output of the A/D is approximated by 
analog integration. This is an excellent approximation since the number 
of samples in 3t/2 is sufficiently large at long ranges. The starting 
time of the range gate integration is designated t $ . 

The signal part of the presum output is designated 

Aj t exp(j 0 j) = S R (L,k) + jSj(L,k) (5) 


where S R and Sj are independent zero-mean Gaussian random variables, 
with variance a g 2 . This is the same value found at the input to the 
presum because of the normalization in our definition. 

The noise components of the presum output are designated 


N R ( L , k ) + 3 Nj(L,k) = 


1 

T 


t +3 t/2 

/ 


t 


s 


N c (t) + j N s (t) dt 


( 6 ) 


where N R and Nj are independent Gaussian random variables with variances 




( 7 ) 
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Formulation of the CFAR threshold for the designated mode in 
search is shown in more detail in Figure 2, In particular, the outputs 
of the DFT doppler filters are given by 


F(m) = |lPSI(k) + j PSQ(k)J exp 
k*0 L “ 


( 8 ) 


for both the early and late range gate outputs. Note that the doppler 
filter outputs from the early and late range gates are correlated; this 
affects the evaluation of the variance of this CFAR threshold. 

The CFAR threshold, T, is formed via the following average 

A ti 

t ■ c i y n p* 1 ’" 1 ’ 1 ) (9) 

• 1=1 L=-l m=0 


where 

I = the RF frequency index, I = 1 I nia)( , I max = 5 

L = the range gate index, L=-l : early gate; L=+l : late gate 

k = the time or pulse index, k = 0,.,., N-l, N = 16. 
m = the frequency of doppler filter index, m = 0,..,, N-l. 

In the next section, the results of the statistical analysis of 
T are discussed. 


t / * * 


F(-l,m) | F{-1 ,m) 



Figure 2. formulation of Shuttle Ku-Band CFAR Threshold for Designated Mode in Search 
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3.0 PERFORMANCE OF CFAR THRESHOLD 

The mean {ensemble average) of the LFM threshold is determined 
in Appendix A and plotted in Figure 3, In this analysis, any doppler 
frequency shift away from the center frequency of the nearest doppler 
filter is neglected. In addition, it Is assumed that the range designa- 
tion is ideal. Both of these assumptions maximize the effect of the sig- 
nal on the statistics cf the CFAR threshold. 

Under the above assumptions, the average CFAR threshold value 
• for the designated mode versus the SNR at the output of the presum is 
plotted in Figure 3, where 

SNR - o s 2 /o n 2 (10) 

As expected, at small values of SNR, the effect of the signal 
from the target disappear and the threshold becomes the value correspond- 
ing to noise only. 

Also included in Figure 3 is the result of the simulation 
reported in [1], The results of our analysis are normalized in Figure 3 
so that the average threshold value at 0 dB coincides with that in [1], 

At this time, we have no explanation for the significant dif- 
ference between the exact analytical and the simulation results since 
their target dependence is less, even though we assumed maximum target 
dependence. 

In Figure 4, the normalized standard deviation of the CFAR 
threshold is plotted versus the SNR at the output of the presum. This 
exact analytical result shows a significant dependence on SNR and a neg- 
ligible dependence on the normalized noise correlation coefficient p. 

For the actual case, we have, as described for the range gates in Fig- * 
ure 1 , ' ' 

P = 2/3. 

Inspection of the plot in Figure 4, however, shows little variation as 
p varies from 0 to 0.7 & 

It is worthy to note the substantial variation in the CFAR 
threshold, particularly at large values of SNR where the normalized > 
standard deviation is greater than 0.3. No attempt has been made to 
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Figure 4. Normalized Standard Deviation of 
Threshold versus Signal -to-Hoise Ratio for 
Various Values of p 
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determine the effects of this variation on the probability of detection. 

On the optimistic side, the threshold setting and the target return are 
correlated; this leads us to conjecture that this variation may not 
appreciably affect the probability of detection. On the pessimistic side, 
there is a substantial variation of the CFAR threshold setting away from 
that developed from the noise-only condition, 


APPENDIX A 


Here we wish to derive the expected value ECT) of the random 
variable T, defined as (see Figures 1 and 2): 


'max 1 N -' 


i I IIU A 

T e V 1 II 

I = 1 L=-l m=0 


0 ) 


where the dependence of F on the frequency range I is explicitly shown 
in (1) and where is a normalizing constant. 

Before we proceed, let us first list the assumptions entailed 
in the following derivations. 

Assumptions 

(1) The in-phase and quadrature-phase components of the sig- 
nal and noise are zero-mean, independent, Gaussian random variables (rv). 

(2) For different RF frequencies, all rv's are independent. 

We can therefore confine our interest to one specific frequency. 

(3) For the same range gate (L) and different time slots (k), 
the noise variables are independent, i.e., 

N^L,^) lN R (L,k 2 ) for k 1 f k 2 

NjO-.kj) J_ Nj(L,k 2 ) 

(4) Real and imaginary parts of either signal or noise are 
always independent, i.e., 

Nr J_ Nj , Sr 1 s n • * 

(5) For the same k but different gate, the noise components 
are correlated. The covariance matrix is 


R =|cov ,N R (-1 ,k),Nj(-l ,k)i,N R (l ,k) ,Nj(l ,k)ji = aj 


1 0 P 0 
0 1 0 p 

P 0 1 0 

0 p 0 1 


A2 

Therefore, for different k's, the noise rv's are independent, regardless 
of the value of L. 

(6) For both range gates and all times k associated with one 
RF frequency, the parts of the signal are Identical, i.e., 

S R (-Uci) * S R (-l,k 2 ) - S R (l,k 3 ) ■ S R (l,k 4 ), 

for every k-j,k 2 ,k 3 ,k 4 [0,N]. Likewise for the imaginary parts of the 
signal Sj. 

As an immediate result of assumption (2), the frequency 
dependence can be dropped and T can be written as 

1 N-l 

T ■ c rW 2 2 i F < L ’ m) i 

L=-l iifO 
N-l 1 

= c Y y | F(L ,m) | 
m s 0 L=-i 
N-l 

- C-l X (in) (2) 

m=o 

where C = C,.r and we have defined the rv X (m ) by 
i max 

X(m) = ! F(-l ,m) | + |F(l,m)| (3) 

In the following, we shall derive the expected value of X(m). From (3), 
it follows that 

E{X(m)} = E{|F(-1 ,m) | } + E{|F(l,m)|} (4) 

■ - <■ 1 ' 

and, because of the symmetry existing between the two gates, (4) simpli- 
fies to 


E{X(m)> = 2E{ |F(1 ,m)|} 


(5) 


It is obvious from (5) and the assumptions made that the first moment 


A3 


(ensemble average) is not a function of the range gate. Henceforth, we 
will drop the gate-index dependence from the symbols listed before in 
order to simply matters. 

Let us call 


y(m) - l F (m ) | 



’.>■•) + F, 2 (m) 


( 6 ) 


so that, from (2), (5) and (6), 


N-l 

E{T) = 20 £ E{y (in) > (7) 

m=0 


Since F(m) is the output of a DFT filter whose input is the set 
{ PSI (k ) + jPSQ(k)}, we have that 


F (in) 


'i- 1 

2 


k=0 


(PSI(k)+jPSQ(k)}e 


• 27rmk 
3 N 


which means that 

N-l 

F R (m) = ^ PSI(k) cos ^-+ PSQ(k) sin (8a) 

k=0 

N-! 

F I (m) = ^ - PSI(k) sin + PSQ(k) cos (8b) 

k=0 v 

I i * s ' t 

Since PSI(k), PSQ(k) are Gaussian rv's, (both signal and noise are Gaus- 
sian), so are their linear combinations F R (in) and Fj(m), which are also 
zero mean. 

We will need the covariance matrix of (F R (m),Fj(m)}. From 
(8a), we have that 



e{f r 2 (i)i)} = & (psi(k) cos £jj»k + PSQ(k) sir/ 


, • „ 2irnik 


fM N-l f 

+ 2 2 ( PSI (k) cos ^Ih. + psq(k) sin pf 

k=0 t=0 ' 
kj^T 

. ^PSI(r) cos + PSQ(t) sin pp^j 


But 


and, similarly, 


E{pSQ 2 (k)| = cr s 2 + a 


A4 


(9) 


E{pSI 2 (k)} - E{s R 2 (k) + N R 2 (k) + 2S R (k) N f (k)} = o g 2 + 0 * (10a) 


(10b) 


For k / x, 

E{PSI (k ) • PSI(t)} 


= E|(S R (k)+N R (k)) (S r (t ) + N r (t ))| ! (Assumption 1) 

= E|s R (k)-S R (t)} + E{N R (k) N r (t)| (Assumptions 3,6) 

(1U) 


* 0 


and, similarly, 


E{ PSQ ( k ) - PSQ(t)} = a 


2 


(lib) 


Also, from Assumption 4, 


E{ PSI (k )- PSQ( k ) } = E{ PSI ( k ) * PSQ (t ) } = E{PSI (t) * PSQ(k)} = 0 


( 12 ) 




Substituting (10), ( 11 ), and (12) into (9) and after some manipulations, 
we get 

e { f R 2 <"')} * (“sH 2 )-? (cos 2 ip/ s m 2 ajfc) 

L'bH « ' 


N-l N-l 


"s 2 I I[f os T < k - T > 


k =0 t = 0 
kj^t 


E { F R 2 (m)} = N(cr s 2 + cr n 2 ) + 2c s 2 ^(N-5) cos (13) 

5=1 

The corresponding result for E\Fj 2 (m) i is easily shown to be the same as 
in (13). ' 

Before we examine (13) closer, let us calculate the cross- 
covariance term: 


l N-l N-l . 

{F R (m)- F I (m)} = E £ (psi(k) cos ^ + PSQ(k) sin 

(k=0 t=0 ' / 

x (-PSI(t) sin ^-+PSQ(t) cos (Assumption 4) 

i ... . . ) 


{ N-l N-l . 

^ ]T (-psi(k) psi(t) si 

k =0 t =0 ' 


» 2wmT rnc 2irmk 

n N C0S “IT" 


+ PSQ(k) PSQ(t) cos sin ^ 


N-l N-l 

.-/n «i" (o-x) 

k =0 t =0 


A6 


We now notice in (14) that, to each term k-r>0 in the double summation, 
there exists a corresponding term -? 3 x-k<0 which, because of the sine 
function being an odd function, cancels with .the former term. Hence, 


E|F R (m)*Fj(m)| = 0 , for all me[0,15] 

and this holds independently of N. The covariance matrix of F R (m), 
Fj(m) can now be written as 


(15) 


A - cov|F R (m),F I (m)| = N 


A(m)*a s 2 + cr n 2 


0 


A(m)<cr s 2 + a n 2 


(16) 


where, from (13) and (16), it follows that 


N-l 


A(«,) = 1 + 2 ^ (l - ft) 


cos 


5=1 


(17) 


We will briefly examine A(m). We have that, for m=0, 


N-l 

HO) = 1 + 2 2 (l - §) = 1 + 2 
5=1 


N-l ' 
(N-l) -I J 5 
5=1 . 


= 1+2 


N-i -%W.n. 


(18) 


To study the case m/0, let us assume that N is divisible by 4. Then, if 
we call = N - £, it follows that , ,,, 


cos 


— = cos pp (IK) = cos 


so that (17) reduces to 
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A(m) = 1 + 2 


Pi/2-1 , 

2 ('- cos2fflt + £2|« 

5=1 ' 


or 


N/2-1 


A(m) = 1+2 cos + cosmir 


5=1 


If we now call 5' = N/2-1, it follows that 


cos 2jTO£l „ cos ZrnmzSl = H) "i CQS 2™i- 


(19) 


so that (19) yields 


A(m) = 1+2 


N/4-1 


v+-i 

£ (l+(-l)' n )cos2« + cosf 

.5=1 


+ cosmir 


( 20 ) 


An immediate conclusion of (20) is that A(m) = 0 for m = odd, For m = even, 
(2) reduces to 


A(m) = 2 


N/4-1 

1 + (- l) 1 ^ 2 + 2 ^ 

5=1 


cos 


2?rm5 


( 21 ) 


The above holds for an arbitrary N divisible by 4. In the general case, 
(21) might yield a nonzero value of A(m), which nevertheless will be 
small. However, it is straightforward to show that, if N is a, power of 2 ^ 
(which is almost always the case in practice), A(m) of (21) identically 
vanishes. Below we summarize the conclusion for the covariance matrix A: 


A8 


A B cov|f r (ui),Fj(iti)} 


No J* + o ^ 0 

® n 

0 cr ^ + 'cr *" 
s n 


v 0 


; ill B 0 


(22) 


i m + 0 


From (22), some useful conclusions can be drawn. First, we notice that 
the noise components at the input affect all doppler filter outputs, 
while the effect of the signal is confined to the m = 0 doppler filter 
only. Furthermore, the effect of the signal on that term is enhanced 
by a factor of N as compared to the noise, Hence, for the moderate-to- 
high signal-to-noise ratio (SNR) environment, we can claim that the 
zeroth filter output is produced by the signal only and the other out- 
puts by the noise only. Finally, we notice that, in all cases and for 
every doppler filter output, the real and imaginary parts of the output 
are independent. This enables us to conclude that y(m) of (6) is a 
Rayleigh rv with mean 


E{y(m) } = Y(ra) 


(23) 


where 


Y(m) =■ 


|J N V + H 


4,7 


; m = 0 


; m i 1 0 


(24) 


From (7), (23) and (24), it follows that 


E{T> = C/2ti- 


*i N V + Ni n 2 + 


or 


(N-D 


E{T> = C-N/27 


c s 2 + °Vn + M " , - 


/N n 


(25) 
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2 2 

or In terns of the sign^l-to-noise ratio SNR ■ cr s /a n » 



r- 


E{T} « C./27.N*cr n 

Ism + i 

t K rfTj 



Typically, the constant in (1) is defined as 

C<| 3 (27) 

In this case, C = 3 1/2N, and (26) modifies to 

I maX 



For the specific application considered, N = 16, so that the factor 

1/N s 0.06 can be neglected for even very moderate SNR (say, SNR>0 dB). 

In this case, we conclude that EtT) varies linearly with (SNR)^. 

0 

A final comment pertains to the values of SNR and cr^ of (28). 
These are the values of the signal-to-noise ratio and noise power at the 
input of the doppler filters. Since an A/D converter precedes these fil- 
ters, the values of these parameters at the input of the A/D converter, 
denoted here by (SNR)^ and a n ^, respectively, relate to the A/D output 
parameters (which are the inputs to the filters) by [2]: 

SNR. 

SNR = -j — - (29a) 

A /° , . 

and 

% = < 29b > 

where 


L A/D 


1.0129 +0.0129 • SNR i 


(30) 








APPENDIX B 


'♦ 2 

We evaluate herein the variance Oy of the* random variable T, 
as defined in equation (1) of Appendix A. 

From (2) of Appendix A, it follows that: 


where 



ril-1 N-1 N-l 

) X(m 2 ) 

m *0 m -|*0 m 2 *0 

m^nig 



( 1 ) 


X(m) =^(-1,111) + F I z (-l,m) + jT R z ( l,m) + F I Z (1 ,m) (2), 

To simplify matters, We will examine two distinct cases: moderate- to-high 
and low SNR. 

1.0 M0DERATE-T0-HI6H SNR 

* 

According to the previous comments of Appendix A, we can justi- 
fiably assume in this case that the zeroth doppler filter output X(0) is 
produced by the signal part of the DFT input only, while all other X(m)> 
mjK) are produced by noise only. 

1.1 Signal Output (m=0) 

E{X 2 (6|} = (from (2))« e|f r 2 (-1,0) + F i 2 (-1,0) + F r 2 (1,0)+F i 2 (1,0) ^ 

/ / 

F r Z (1.0) + F, Z (1, 

Since we have assumed that all random variables appearing in (3) are 
exclusively signal functions, we can further invoke Assumption 6 of Appen- 
dix A, according to which W 

*5$ 

and Fj(-1 ,0) = Fj(l ,0), © ' 

m 

© 

a 


F R (-1,°) = f r (i,o) 


i I'M 





B2 




Substituting the above Into (3) yields 

?{x 2 (0)} ■ 4e{f r 2 (-1,0)} + 4e{Fj 2 (-1,0)} (4) 

2 2 

From (22) of Appendix A, we have (setting a » a) that 

5 n 

e|F r 2 (-1 ,0)j • e{f, 2 (-1,0)} - N 2 o s 2 (5) 

so that (4) and (5) combine to give 

e{x 2 (0)} * 8N 2 a s 2 (6) 


1.2 Noise Outputs (n^O) 

E Uo ! ” ! ’ " E ( r R 2( " 1 '” l) + F I 2( * 1 ’ n,) + f r 2 ( 1 - 1 ") + F i 2 0.m) 


+ 2j( F R 2 (-l .in) + F, 2 (-l .ni)) ( Fr 2 (1 ,m) + F, 2 (l ,m)jj (7) 


Since all random variables are produced by noise, it follows from (22) 
of Appendix A that 

e{f r 2 (-1,«,)} » e{Fi 2 0 ,m)} ■ e|f r 2 (1 ,m)j = e{f, 2 (1 ,m)} « N^ 2 h o 2 

so that 

e{ X 2 (m)l = 4Na 2 + 2 D (8) 

ntfO ’ n 

where < • * " 

D - e'|JP r Z (-1 ,m) + F I 2 (-l,iTt))(F R 2 (l,i») + Fj 2 (1 ,m)) | (9) 

We now proceed to evaluate D. To <IP‘this, we need the joint statistics 
of the zero-mean, Gaussian rv's appearing in (9). Let us find the covar- 
iance matrix Rof {F R (-l,m), Fj(-l,m), F R (l,m), Fj(l,m)}, m/0. 


D3 


From (8a) of Appendix A, we have that (assuming only noise Input), 
, s (p-1 

E\F^(-l,m)? * E< (n r (k,-1) cos + Nj(k,-1) sin 

* ‘ (|k*u V 


(because of assumption 4) 

N-l 

* a n 2 2 ( cos2 nin +s, * n2 hp ) * >% 2 * ° 2 no ) 

I A * 


and similarly for the others. 
Furthermore, 

< > fc 

EjF R (-l,m).F R 0,n>)} - E{2, Z 
^ ' (k*o t=o 


! \ iiri n-i . 

F r (-1 >m)'F R (l ,m)> > £ (N R (-l,k) cos Nj(-1 ,k) sin 

' (k*o t«o ^ 

(n r O,t) cosSp + Njd.T) slnSf-)' 01) 


which, because of Assumptions 4 and 5, reduces to 


e!f r (-1 ,m).F R 0 ,m)j = e£ N R (-l,k)N R (l,k) cos 2 

* ‘ IkcQ 


+N T (-l,k)*N T (l,k) sin 2 ^"4 * N-« 2 -p “ pd 2 (12) 


Hence, the covariance matrix Ris 


R * cov<F R (-l,m),F I (-l,m),F R (l,m),F I (l-m)> = d 


I 4 

_ 1 


p 0 1 0 

0 P 0 1 


1 0 p 

0 1 0 p ( 13 ) 


Comparing (13) and Assumption 5, we draw the interesting conclusion that 
the noise statistics remain unaltered after passing through the complex 
DFT doppler filters. 
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We are now in a position to evaluate D. If we change to 
polar coordinates 

F R (-l,m) = V £ cos i> E F R (l,m) = v, cos 

and 

= V E sin * E Fj(l ,m) = v L sin * L 

then (9) is simply written as 


04 ) 


D = E 




05) 


The jnlnt distribution of the envelopes II, and v L . for the four Gaussian 
rv s coi related as in (13), is found from [3] to be 



06 ) 


where y.) is the zeroth-order modified Bessel function. From (15) and 
(16), it follows that 


■// 


v 2 v 2 

V E V L 


D * jj -pf ? - r oi 

0D K 


ToTrj- «p 


u 2 (V E 2 tV L 2 ) 


2|R 


T7|— ? dV £ dV L < 17 ) 


If we make the transformation 


2 | R | ^/ 2 

V E = 2 — *yexp {2z> 


2 |r i V2 

"l = — • y exp {-2z} 

a 
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whose Jacobian is 


and, upon substituting into (17), we get 


i ,r i r 

VI 


dy • y d • i n (yp) 


■(£*• 


exp{-y cosh (2z )}) (18) 


The inner integral equals Kg (y ) [4], where K Q (y) is the zeroth-order 
modified Hankel function. Therefore, from (18), 


/ y 2 J Q (jpy) k 0 (y) dy 

° J o' 


where we have used the fact that I Q (n) = J Q ( jn ) (J q( • ) is the zeroth-order 
Bessel function). 

The integral of (19) can be evaluated. From [5], we find that , 


j(az) J g (bz)-z“ Y dz 




, 20 ) ; 


where £F.j(*) 1S ^e Gaussian hypergeometric function [5]. Applying (20) 
for a = 0, 8 = 0, y= -2, a = 1, b = jp,we get from (19) that 

D = i_Ri . | . 2 F 7 (3/2; 3/2;l; P 2 ) (21) 
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From [5], we find that 


2 F,(3/2s 3/2; 1; p 4 ) - (l-p 2 )‘ 2 - 2 F,(-l/2; -1/2; 1; p 2 ) 

= 0-p 2 )' 2 -|e(p 2 ) - |(l-P 2 )-k(p 2 )] (22) 


where k( * ) and E(-) are the complete elliptic integrals of the first and 
( second kind [6], 

Upon substitution of (22) into (21), we find that 


D = 


T 


(1-p*) 


-2 


2E(p 2 ) - (1-p 2 ) k(p 2 ) 


(23) 


From (13), the determinant | R | can be evaluated to be 
| R | = c?(l-p 2 ) = N 4 P n 8 (l-p 2 ) 


(24) 


From (23) and (24), we finally get 


D = 



2E(p 2 ) - (1-p 2 ) k( P 2 ) 

J 


(25) 


From (8) and (25), we derive that 





2 + 2E(p 2 ) - (1-p 2 ) k(p 2 ) 


(26) 


Equations (6) and (26) yield the result for the first summation in (1). 

The second term (double summation) is now evaluated. 

First let us assume that m^ ^nig and none of them equals zero. 

In this case, both X (m.|) and X (mg) are produced by noise, and from (2), 

'it 




EjXOn,) X (m 2 )> = fy(F R c (-l ^ ) + Fj fc (-1 ^ )) .(F R fc (-l ,m 2 ) + F^-l ,m 2 )tf 

+ 4/(^0 ' m i ) + Fj 4 (1 , 111 , )) -(f r 2 (1 ,m 2 ) + Fj 2 (1 ,ntj) ) ) 

+ e|J(f r 2 (-1 .m, ) + F, 2 (-l ,m, )) .(f r z 0 ,m 2 ) + F, 2 (l ,m 2 ))| 

+ f[/( F r 2 (1 ,iii, ) + F[ 2 (l ,in, )) (f r 2 (-1 ,m 2 ) + Fj 2 (-1 ,m 2 ))| (27) 

We will now show that all random variables under the same square root are 
mutually independent. To do this, let us evaluate 

( ) It? t* , 2Trnii k 2irm 1 k\ 

EjF R (-l,m 1 )-F R (-l,m 2 )> = eO ^ ( N R (-1 ,k) cos -tf-* Nj(-1 ,k) sin 

\ m a!m ) / 1<=0 T = 0 


m^m 2 


( 2irmpT 2irmpX 

N r (-1,t) cos — pj — + Nj(-1 ,t) sin — j- 
(Assumption 4) 

\<r-> n 2irmi k 2irm ? k 

= eO N r (-l,k) cos — cos — fl— 


+ Nj (-l,k) sin 


2irm,k 2TT|Tlpk 

! *• *» irk _ 


■ % 2 I 

k=0 


2irk(m^ - m 2 ) 


For N a power of 2 and m^m 2> it is easy to show that the summation in 
(28) vanishes which, in, turn, means that F R (-l,mj) and F R (-l,m 2 ) are inde- 
pendent. Identically, we show that Fj (-1 ,m-j ) J_F j (-1 ,m 2 ) . Similarly, 
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N-l 


E K (- l » n , 1 )* F R 0 .m2)} "% 2 S cos R 

k=0 


2nk(ni 1 -m 2 ) 


= 0 , 


Finally, by interchanging 1 and -1, we get the symmetrical results for 
both gates, which justify our earlier claim that all rv's in (27) are 
independent. Therefore, 


E < X (m-| ) X (nig ) 
^ O/m-jAitg/O ' 


4E 2 j Jf^(-1 ,'«! ) + Fj 2 (-1 , m i ) j 

4 -(*M - 


(29) 


Second, assume that m 1 = m 2 =m^0. Because X (0) is produced by the sig- 
nal, X (m) by the noise and, due to the high SNR assumption, we get 


E{ X (m) X (0)} = E{ X (m)} -E{ X (0)} (see (22) of Appendix A) 
m/0 

= fU *jN 2 a s 2 • /2? -JnaJ = 2 N/N 0 ^ (30) 

From the N 2 terms involved in (1), one term is given by (6), (N-l) terms 
are given by (26), 2 (N-l ) terms are given by (30) and (M-l ) * (N-2 ) terms 
are given by (29). Summarizing, we get that 




8N 2 a s 2 + (N-l)-2Ncr n 2 [2 + 2E(p 2 ) - (1-p 2 ) k(p 2 j 
+2(N-1)-2 ttN/N a s <r n + (N-l ) (N-2) 2 -itN^ 2 


2C 2 N 2 cj 2 
n 


4.SNR+(l-l) [2 + 2E(p 2 ) (1-p 2 ) K(p 2 ) 


MMlj SNR + (i 1) . (N-2) •ir 
/N v 


t / * * 


+ 


(31) 
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But from (28) of Appendix A, we have that (for high SNR) 


E{T} = C • N/cJtF • 


/SNR + 


( M ) 

\ At) 


therefore , 


E 2 {T} 


2C 2 N 2 • 


, SNR + x(h& - + 


(32) 


From (31) and (32), it follows that 


var{T) = a. 


= E{T 2 } - E 2 {T) 


= 2C 2 N 2 a n 2 


(4-ir) SNR + (l J 2 + 2E(p 2 ) - (1 -p 2 ) K(p 2 ) 


+ (N-2)ir - (N-I)tt 


or 


2 2 2 
oj = 2(TN 


a„ 2 [(4-ir) SNR + (l -J-) H(p )] 


where 


H(p) = 2 - it + 2E(p 2 ) - (1-p 2 ) K(p 2 ) 


(Moderate-to-high SNR) 


(33) 


The primary quantity of interest is the ratio cr-j-/E{T> which, 
for case 1.1 is given by 


° T _ gC» <yJ(4-») SNR * (l-^Htpj 
E{T} CN n /& 'G n (/SNR + *tl\ 

0 /N / 



BIO 



(Moderate-to-high SNR) 


VERY LOW SNR (OR SIGNAL ABSENT) 


In this case, all the doppler filter outputs are produced by 
noise. If signal is also present (in order for this assumption to be 
valid, even for the m = 0 filter), we must have that 


Ncf « o c « 


For N=16, such an inequality is satisfied for SNR in the order of 
-20 dB or less. 

From (26) (which also holds for m = 0 now) and (29) (which is 
now restricted only to m^m 2 )» we get that 


e[t 2 } = C 2 ■ 2N 2 o n 2 fa + 2E(p 2 ) - (1-p 2 ) k(p 2 )]+ N(N-l) 2* n 
■ 2C 2 N 2 o n 2 [2 + (N-l) + 2E(p 2 ) - (1-p 2 ) k(p 2 )] 

For very low SNR, we have from (28) of Appendix A that 


E{T> = C • . N • 


E 2 (T} = 2 it C 2 N 3 a n 2 


therefore 



and 
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var{T} = a T 2 = 2C 2 N 2 a n 2 2 + (N-l )ir + 2E(p 2 ) - (1-p 2 ) k(p 2 ) - Nir 


or 


var{T} = 2C 2 N 2 o n 2 • H(p) 


(36) 


, where H(p) has been defined in (33). The ratio cr^./E{T> for this case 
(2.0) is found to be 


ETTT 


^■NCo n 4TTpy 
/l -C-N4T.a n 


°T _ 

/H(p) 

E{T) 

si Ntt 


(Very low SNR) 


(37) 


independent of the SNR, as expected, The ratios of (36) and (37) have 
been plotted in Figure 4 for various values of p. The function H(p) 
depends on the complete elliptic integrals E(p ) and K(p ), which have 
been tabulated for different values of the argument p 2 [6]. In Table 1, 
we give some values of H(p) as a function of p. 


* * * ♦ 
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